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A realistic five-dimensional warped scenario with all standard model fields propagating in the bulk 
is proposed. Mass hierarchies would in principle be accounted for by judicious choices of the bulk 
mass parameters, while fermion mixing angles are restricted by a A(27) flavor symmetry broken on 
the branes by flavon fields.The latter gives stringent predictions for the neutrino mixing parameters, 
and the Dirac CP violation phase, all described in terms of only two independent parameters at 
leading order. The scheme also gives an adequate CKM fit and should be testable within upcoming 
oscillation experiments. 
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I. INTRODUCTION 


The understanding of flavor constitutes one of the most stubborn open challenges in particle physics [l|. Two 
aspects of the problem are the understanding of fermion mass hierarchies as well as mixing parameters. Various types 
of flavor symmetries have been invoked in this context . These efforts have been partly motivated by the original 
success of the tri-bimaximal mixing ansatz [ir)| . The resulting non-Abelian flavor symmetries are typically broken 
spontaneously down to two different residual subgroups in the neutrino and the charged lepton sectors, lea ding to 


zero reactor mixing parameter, 613 = 0 . However, the measurement of a non-zero value for the reactor angle lll-ll4 


implies the need to revamp the original flavor symmetry-based a ppr oaches in order to generate 613 ^ 0 [l5| or else 


look for alternative possibilities, such as bi-large neutrino mixing 16l-ll8l|. 




The existence of warped extra-dimensions has been advocated by Randall & Sundrum 19| as a way to address the 
hierarchy problem, since the fundamental scale of gravity is exponentially reduced from the Planck mass down to the 
TeV scale as a result of having the Higgs sector localized near the boundary of the extra dimensions. Moreover, if 
standard model fermions are allowed to propagate in the bulk and also become localized towards either brane, the 
scenario can also address the flavor problem possibly acting in synergy with the flavor group predictions. This is what 
we do in the present paper. 

The idea of combining discrete flavor symmetries and extra dimensions is quite attractive and has already been 


discussed in the literature within the context of large extra dimensions 20l-l22l|. warped extra dimensions 23l-l27| and 






holographic composite Higgs models 28|-l30|. However, such models try to generate tri-bimaximal neutrino mixing. 


which ha s been ruled out by the measurement of the reactor angle ^13 
data 
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I4| and also global fits of neutrino oscillation 


3l|. One of us has constructed a warped extra dimension model with S 4 flavor symmetry where democratic 




mixing is produced at leading order and non-zero ^13 can arise from subleading corrections [32| . In this work, we shall 
re-consider the issue of predicting flavor properties in particle physics by combining the conventional predictive power 
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inherent in the use of non-Abelian flavor symmetries with the presence of warped extra-dimensions. We propose a 
warped five-dimensional scenario in which all matter fields propagate in the bulk and neutrinos are treated as Dirac 
particles. Our model can accommodate all the strengths of the standard model Yukawa couplings and resulting fermion 
mass hierarchies by making adequate choices of fermion bulk mass parameters, while the fermion mixing parameters 
can be restricted by means of the assumed flavor symmetry. We present a A(27) based flavor symmetry which nicely 
describes the neutrino oscillation parameters in terms of just two independent parameters, leading to interesting 
correlations involving the neutrino mass hierarchy and the leptonic Dirac CP phase, not yet reliably determined by 
current global oscillation fit Our predictions include a neat leading order relation between the solar and reactor 
mixing parameters which should be tested at future oscillation experiments. 


II. BASIC STRUCTURE OF THE MODEL 


In this section we present the basic setup of a warped five-dimensional (5D) model for fermions, constructed under 
a A(27) 0 Z 4 0 Z 4 flavor symmetry. The 5D field theory is defined on a slice of AdSs, where the bulk geometry is 
described by the metric 


ds^ = e — dy'^ , 


( 1 ) 


with = diag(l, —1, —1, —1) and k as the AdSs curvature scale. The fifth dimension y is compactified on SijZ 2 , 
and two flat 3-branes of opposite tension are attached to the orbifold fixed points, located at y = 0 (UV brane) and 
y = L (IR brane). 

The electroweak symmetry of the model is promoted to Gbuik = SU{2)l 0 SU{2 )r 0 U{1)b-l in order to avoid 
excessive contributions to the Peskin-Takeuchi T parameter 33, |^. The gauge group Gbuik breaks down to the 
standard model electroweak (EW) group Gsm = SU{2)l®U{1)y on the UV brane by the boundary conditions (BCs) 
of the gauge bosons. Furthermore, a bulk Higgs field with (5't/(2)b, SU{2)r) quantum numbers 


iJ- (2,2) 


( 2 ) 


is responsible for the spontaneous symmetry breaking (SSB) of Gsm- The 5D Higgs field H{x^, y) can be decomposed 
into Kaluza-Klein (KK) modes as 


H{x^ ^y) = H{x^) heavy KK Modes. 

yL 


For an adequate choice of BCs, its zero mode profile fnijj) can be written as 


( 3 ) 




' 2fcL(l-/3) 

1 _ g-2(l-/3)feL 


gfcLg(2-/3)fc(y-L) 


( 4 ) 


where we have introduced the Higgs localization parameter /3 = -^4-1- in terms of the Higgs field bulk mass 

parameter niH- In the present work, we assume that the vacuum expectation value (VEV) of the Higgs zero mode is 
of the form 

<"<"'>>-5(0 O' 

and it is peaked toward the IR brane, allowing for a TeV scale EW SSB and inducing the Gbuik breakdown to 
SU{2 )d 0 U{1)b-l on that brane. 

Three families of fermion fields are required to describe each generation (labeled by z = 1,2,3) of quarks and 
leptons. All fermion fields propagate into the bulk and transform under the minimal representation of the gauge 









3 


group SU{2)l ® SU{2)fi [33|,[3^. In the lepton sector the three multiplets of the model are given as 
/ „[++] \ / „[+-] \ / 1 \ 


= 


J++] 


( 2 , 1 ), 


= 


( 1 , 2 ), 




d+-l 


( 1 , 2 ), 


( 6 ) 


while for the quark sector we have 

d 


^Qi = 


j[++] 


( 2 , 1 ), 


,[+-] 




( 1 , 2 ), 




^[+-1 


( 1 , 2 ). 


(7) 


Notice that we have a separate SU{2)ii doublet for every right handed fermion. In the above equations, fields with 
different sign assignments must be understood as independent. The bracketed signs indicate Neumann (+) or Dirichlet 
(—) BCs for the left-handed component of the corresponding held, on both UV and IR branes. The right-handed part 
of the held satishes opposite BCs. Only helds with [-1—k] BCs have left-handed zero modes, whereas right-handed zero 
modes exist solely for helds with [-] BCs. The KK decomposition for such helds has the form 


p2ky ^ >, 

, y) = + heavy KK modesj , 




( 8 ) 


-k heavy KK modesj 


with ■)/, = and zero mode prohles 


36-3 


ff\y,CL) = 


{l-2cL)kL 

o(l —2c2.)fcL _ 2^ ^ 


Chky 


/!?’<,/.«)= (9) 


where cl and c/j are the bulk mass parameters of the 5D fermion helds in units of the AdSs curvature k. Thus, the 
low energy spectrum contains left-handed doublets iiL = (r'iL,eiL), QiL = {uiLjdiL), alongside right-handed singlets 
k'iR, em, UiR, diR. In the following, we identify all standard model helds with this set of zero modes (he. the so called 
zero mode approximation, ZMA). For future convenience, we denote the havor components of charged leptons and 
quarks as 61 ^ 2,3 = e, /r, r; Qi, 2,3 = U,C,T; mi_ 2,3 = u, c, t; di^ 2,3 = d, s, b. 

In the present work, we choose the havor symmetry to be A(27), augmented by the auxiliary symmetry Z 4 ® Z'^. 


The group A(27) was originally prop osed to explain the fermion masses and havor mixing in Refs. 39|, l40|, and has 


been used for Dirac neutrinos in 


4l| by one of us. Here we study its implementation in a warped extra dimensional 
theory. The havor symmetry A(27) 0 Z 4 0 Z 4 is broken by brane localized havons, transforming as singlets under 
Gbuik- We introduce a set of havons 5, cri, a 2 localized on the IR brane, and a havon ip localized on the UV 
brane. Both ^ and p are assigned to the three dimensional representation 3 of A(27), while CTi and (T 2 transform as 
inequivalent one dimensional representations lo,i and lo,o respectively. A summary of the A(27) group properties 
and its representations can be found in Appendix There are two different scenarios for the model, determined by 
the two possible VEV alignments for namely: 


(0 = (0,1,0)U5, 

(0 = (l,w,l)?;j. 


Case I, 
Case II, 


( 10 ) 


with u) = As indicated above, we will denote the models described by each alignment as cases I and II 

respectively. Note that the case II vacuum pattern frequently appears in the context of geometrical CP violation 
43|. The VEVs for the remaining havon helds are 


42 


{^) = ( 1 , 1 , 1 )’ 


V ’ 


(c^l) = 


( 0 ' 2 ) = ' 


( 11 ) 


Further details regarding this vacuum conhguration are offered in Appendix |B] 
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III. LEPTON SECTOR 

Once the basic framework has been laid out, we are in position to discuss the structure of the lepton sector and its 
phenomenological implications. As we will show below, charged lepton as well as Dirac neutrino masses are generated 
at leading order (LO), and non-zero values for the “reactor angle” @13 arise naturally. The model is predictive, in 
the sense that the three mixing angles and the Dirac CP phase will ultimately be determined in terms of only two 
parameters. 


A. Lepton masses and mixing 


Field 



'I'm 

■i’r 


4'.2 

4- ,3 

H 



CTl 

(T2 

A(27) 

3 

lo,o 

ll,0 

l2,0 

lo,o 

lo,o 

lo,o 

lo,o 

3 

3 

lo,l 

lo,o 

^4 

1 

1 

1 

1 

-1 

i 

-1 

1 

1 

-1 

1 

i 


1 

i 

i 

i 

-1 

-1 

-1 

1 

—i 

1 

-1 

-1 


TABLE I: Particle content and transformation properties of the lepton and scalar sectors under the flavor symmetry A(27) ® 

^4 » z'^. 


The transformation properties of leptons and scalars under the family symmetry A(27) (D Z 4 ® Z 4 are given in 
Table m Note that the Higgs field is inert under the flavor symmetry. Since the three left-handed lepton doublets are 
unified into a faithful triplet 3 of A(27), they will share one common bulk mass parameter ci- On the other hand, both 
right-handed charged leptons and right-handed neutrinos are assigned to singlet representations of A(27). Therefore, 
there are six different bulk mass parameters Ce^ and Cv^ (i = 1,2, 3) for these helds. From the particle transformation 
properties we can write the most general lepton Yukawa interactions that are both gauge and flavor invariant at LO: 


Vg f _ ~ ~ _~ 

+ I 2/11 + 2/31 + 2/22 


+2/i3($cri4'^) +y 33 {^a*-$e) \s{y- L) + h.c. 


( 12 ) 


with H = T 2 H*T 2 , and Ti as the Pauli matrices. After electroweak and flavor spontaneous symmetry breaking, all 
leptons develop masses dictated by the above Yukawa interactions. The generated masses are modulated by the 
overlap of the relevant zero mode fermion profiles, the VEV profile of the Higgs, and the flavon VEVs given in 

Eqs. (nnmn). 

From Eq. (Ha, The mass matrix m; for charged leptons is 


mi = 


' ” -rm 


{LK)i A ^2 

where Ui stands for the so-called magic matrix 





/ill 

1 UJ 
\ 1 


(13) 


LO 


( 14 ) 
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and ye,^,T are modified Yukawa couplings defined as 

ye,p,r — ^ei) ) 

in terms of the overlapping function 

F{y,CL,CR) = fL\y,CL)fR\y,CR)fH{y) 


(15) 


2(l-/3g)(l-2cL)(l + 2cfl) 

1^]^ _ g —2(1—/3H)feij 'yg{l — 2cL)kL _ |^g(l+2cK)fcL _ 


^-{l-0H)kL (2-PH-CL+CR)ky 


(16) 


Given that Cf/C/j = 1, the diagonalization of the charged lepton mass matrix is straightforward, leading to charged 
lepton masses of the form 




ye,^,T y 


(17) 


(LA)i AV 2 ' 

Analogously, taking into account the two distinct VEV alignments for the flavon triplet ^ in Eq. (ITUl) . the neutrino 


mass matrix for each respective case can be written as 


T 1 Uf i; 

m„ = 


{LA')I A' P2 


I yiii^ 0 2/131^ 

0 ^ y 22 ^ 0 ^ 

\ 2/31 0 y33-^ 


II 1 u 

= 


(LA')^ A'y 2 


with 


and 


I yii-K^ 0 yis-^h- 

V3Eo 0^ y22^ 0^ 

Vyaiy^ 0 y33-^ 

yij ~ yijFi^Ly c^, ), 

/ w 1 1 

A 




1 w 1 

y 1 1 w 

Thus, the diagonalizing matrix for the neutrino sector can be parameterized as 


/ 


U.= 


cos 9 1 , 

0 


0 sind^e*'^'' 
1 0 


y—sin 0^6 0 cos 9^ 


( cos 01 / 0 s\\3 9ye^'^‘' 


U}} = En 


0 


1 


0 


y—sin 0^6 0 cos 0 i/ 


In terms of the auxiliary functions 

xt = \ y 3 i ? + \ y 33? ± lyiiP ± I013P, E. = ynyss - y * i 3 ^i , = 0112/31+^13033 ^ 

the relevant parameters of the model, 9i, and are given by 

tan 20 j/ = 2 |Z^|/Y“, = arg(u^^Zj/) , 

and the neutrino mass eigenvalues for both NH and IH are determined as 


(18) 


(19) 


( 20 ) 


( 21 ) 


( 22 ) 

(23) 

(24) 

(25) 
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• Case I 


NH: 

mi = ^M (XF,Y„), m2 = V2\y22\, m3 =-^M+(X+, Fj,), for cos 20„ > 0 , 

(26) 

IH: 

mi = ■^M+(X+,Fi,) , m2=U2|022|, m3 = (X+, F^), for X“ cos 20^, < 0 , 

(27) 

• Case H 



NH: mi 

= Y^uiA/“(X+, Fj,) , m2 = V3v2 y22 , m3 = ,Y„), for X“cos20i, > 0, 

(28) 

IH: mi 

= {XF,Yi,) , m2 = a/ 3 u 2 1^22! , m3 = y^viM“(X+,Fi,), for X“ cos 20^ < 0 , 

(29) 

where we have defined 



= \Jx± \fx^ -4|2/|2, (30) 


and 


Va 


1 Va^ V 

{LA')i A' A' V2 


a = 1,2. 


(31) 


Without loss of generality, the angle 9^, is restricted to the interval [0,7r]. Notice that X“cos20,^ = 
2|Zj^| cos^ 20;^/sin2di,. As a result, for non-vanishing values of Z^, the neutrino mass spectrum displays Normal 
Hierarchy (NH) provided 0 < 0^ < 7r/2, whereas Inverted Hierarchy (IH) is realized for 7r/2 < < tt. The angle (p,^, 

on the other hand, can take any value in the interval [0, 27r]. 

At leading order, the lepton mixing matrix t/pMNS = U^Ui, becomes 


r/i _ J_ 

'^PMNS — yg 


cos 9^ — Loe 


ycos^i, —w^e *'^‘'sin0y lo cos sin ( 


(32) 


uF 


PMNS 


— ILU 

75 


^ cosdj/— sindi. 

w cos 01, — sin 0^, 

y wcosdi/— e"®"^" sindj. 


1 

1 


cos 01, -I- e®''®" sin 0^ 
cos 0„ + we®'^'' sin 0„ 
cos 01, + we®'^'' sin 0^, 


(33) 


In both cases, the solar, atmospheric and reactor angles can be written in terms of 0i, and Lp^ as 


sin^ 012 
sin^ 023 
sin^ 013 


1 

2 — sin 201, cos (pi, ’ 

1 — sin20,y sin(7r/6 — p^,) 
2 — sin 201, cos p„ 

i (1 -I-sin 201, cos (^i,) . 

O 


(34) 


A convenient description for the CP violating phase in this sector is the Jarlskog invariant Jqp = 
Im[C/*i U^sU^iUes] |44| . which in this parameterization takes the compact form 


d^CP 


— -j= cos 201,. 
6v 3 


(35) 


It is worthy of attention the independence of Jcp upon and the simple predicted relation between the solar and 
reactor angles 0i2 and 0i3: 


sin^ 012 cos^ 013 = 


1 

3 ■ 


(36) 
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NH IH 




FIG. 1; la, 2a and 3cr ranges of sin^ 612 (green), sin^ 623 (blue) and sin^ @13 (red) for normal (left panel) and inverted (right 
panel) neutrino mass hierarchies. Best-fit contours for sin^ 613 (sin^ 023) are indicated by dotted (short-dashed) lines. The 
long-dashed contour in the left panel represents the local minimum in the hrst octant of 023 . 


B. Phenomenological implications 


As shown above, only two parameters are required to generate the three angles and the Dirac CP violating phase 
characterizing the lepton mixing matrix, making this model highly predictive. In the remaining part of this section 
we explore in detail the predictions for the lepton mixing parameters and the neutrino mass spectrum. 

In Figure [l] the 0^ - ipi, parameter region compatible with experimental data is delimited using the global fit 
of neutrino oscillations given in ^ for each mass ordering, shown as the left and right hand panel. The model 
can reproduce successfully the best-fit values for the atmospheric and reactor angles, reaching simultaneously the 2 a 
region for the solar angle. The intersecting points of the “central” or best fit curve in the sin^ 0i3 contour and the 
corresponding ones in the sin^ 023 contour are located at 


NHi : 

QvIt: 

= 0.204(0.296), 


= 0.924, 

NH 2 : 

9vh 

= 0.193(0.307), 


= 1.031, 

IH ; 

9vh 

= 0.707(0.793), 


= 1.917, 


(37) 


where NHi denotes the best-fit contour of sin^ 023 , and NH 2 corresponds to its local minimum in the first octant. 
Notice that the numbers in parenthesis denote the intersection values within the range 9^, € [7r/4,7r/2] U [37r/4 , tt]. 

Once we have determined On and Lp^ from the central values of the atmospheric and reactor oscillation global fits, 
the predictions for the solar angle and the Jarlskog invariant can be straightforwardly obtained using Eqs. (|M1 [35|). 
For completeness, in Table HIl we present the full set of mixing parameters derived from the points defined in Eq. (|37)) . 

Remarkably, the central prediction for sin^ 0i2 falls very close to its Icr boundary. In addition, notice that the Icr 
range of Jcp is entirely contained in the region 9^, € [0,7r/4] U [37r/4, tt]. 

We conclude this section bringing forth a consistent realization of lepton masses and mixing angles. In the numerical 
analysis, we assume that the fundamental 5D scale is ~ A ~ Mpi, with Mp\ ~ 2.44 x 10^® GeV as the reduced 
Planck mass. We also set the scale K' ~ k' = ke~^^ ~ 1.5 TeV in order to account for the hierarchy between the 
Planck and the electroweak scales, allowing for the lowest KK gauge boson resonances (with masses rriKK = 3 ~ 4 
TeV) to be within reach of the LHC experiments. The Higgs VEV is identified with its standard model value v ~ 246 
















NHi 

NH 2 

IH 

sin^ 023/10 ^ 

5.67 

4.73 

5.73 

sin^ 013/10“^ 

2.26 

2.26 

2.29 

sin^ 012/10-1 

3.41 

3.41 

3.41 

Jop/lO”^ 

-(-b)2.71 

-(-f)3.37 

+(-)2.57 


TABLE II: Central predictions for sin^ 612 and Jcp obtained from the central values of the atmospheric and reactor angles 
reported in Ref. [s^. The sign of Jcp in the parentheses corresponds to the bracketed prediction for 6^ in Eq. pTll. 


GeV, and the ratios v^/A', A', A', are all fixed to 0.1 (thus considering real-valued flavon VEVs). The 

Higgs localization parameter /?, common to all mass matrix elements, is chosen as 0.95 in the following discussion. 

As an as illustrative example, we can choose = 1.85, Ce = —0.27, = —0.44, Cr = —0.71, \ye\ = 0.861, 

\y^\ = 0.898, \yr \ = 0.994 to generate the charged lepton masses me = 0.511 MeV, = 105.7MeV, mr = 1.777GeV. 
For the neutrino sector, benchmark points (BPs) in parameter space are given in Table Hill There, the four BPs are 
labeled according to their hierarchy scheme and case as NH-I, NH-II, IH-I, IH-II. One sees that, indeed, the large 
disparity between charged lepton masses is reproduced for Yukawa couplings of the same order of magnitude. 



NH-I 

NH-H 

IH-I 

IH-H 

C-ui 

-1.40 

-1.41 

-1.39 

-1.40 

Cu2 

-1.38 

-1.40 

-1.33 

-1.35 


-1.34 

-1.36 

-1.34 

-1.36 

011 

-1.000 - 0.307i 

0.282 + 1.166i 

0.752 -b 0.096* 

-0.674 + 0.520* 

013 

-0.451 -b 0.631i 

0.031 - 0.880i 

0.919 - 0.432* 

1.026 - 0.542i 

022 

0.860 -b 0.353i 

0.097 - 1.088i 

-0.905 - 0.194* 

0.974 + 0.431i 

031 

0.667 -b 0.397i 

0.001 - 0.881* 

0.941 -b 0.383* 

-1.070 + 0.450* 

033 

0.792 - 0.683i 

-0.324 + 1.154i 

0.746 - 0.136* 

0.829 - 0.191i 


TABLE III: Benchmark points for the neutrino sector, featuring both NH and IH in Cases I and 11. 


The neutrino masses, splittings and mixing angles associated to each BP are displayed in Table ITVl All the obtained 
neutrino oscillation parameters are consistent with the global fit in Ref. 31|. In particular, the reproduced atmospheric 
and reactor angles lie comfortably in their respective Icr region, whereas the solar angle values are contained in the 
2(7 range, very close to the la boundary. 



NH-I 

NH-II 

IH-I 

IH-H 

mi [eV] 

1.80 X 10-® 

2.59 X 10-® 

4.88 X 10-2 

4.89 X 10-2 

m 2 [eV] 

8.90 X 10-^ 

9.10 X 10-® 

4.96 X 10-2 

4.97 X 10-2 

m 3 [eV] 

4.98 X 10-2 

4.99 X 10-2 

2.41 X 10-® 

3.50 X 10-^ 

Amii [10-®eV2] 

7.60 

7.60 

7.50 

7.48 

|Amii|[10-®eV=l 

2.48 

2.48 

2.38 

2.38 

sin^ 012 / 10 “^ 

3.41 

3.41 

3.41 

3.41 

sin^ 023/10“^ 

5.67 

5.67 

5.73 

5.73 

sin^ 013/10-2 

2.26 

2.26 

2.29 

2.29 

Jcp/10-2 

-2.71 

-2.71 

-2.58 

-2.57 


TABLE IV: Neutrino masses and oscillation parameters associated to the four chosen benchmark points. 
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IV. QUARK SECTOR 


Field 

4'c/ 

'I'c 




4't 

'kd 


^6 

A(27) 

lo ,2 

lo,i 

lo,o 

lo ,2 

lo,o 

lo ,2 

lo,i 

lo,o 

lo,l 

^4 

—i 

—i 

—i 

1 

1 

—i 

1 

—i 

—i 

^4 

1 

1 

1 

-1 

-1 

-1 

-1 

-1 

-1 


TABLE V: Particle content and transformation properties of the quark sector under the flavor symmetry A(27) ^ (gi Z'^. 


The quark transformation properties under the family group A(27) ® ® Z'^ are given in Table |Vl At leading 

order, the most general invariant Yukawa interactions can be written as 

^^G ' 


= 


(A')t 


yuucr2'^uH'i/u + yctCTi'i'cH'i’t + ytc'y2^TH'^c + 
+2/dsCT^^'c/i^^'s + ydbcri'^uH^b + ysdcr^'^cH^ d 
+yssCfi'^cH'^s + ybb(Ji'^TH'l!}}^5{y - A) + h.c. 


(38) 


Again, after spontaneous electroweak and flavor symmetry breaking, the mass matrices for the up and down quark 
sectors read 

( 


1 


(LA')2 72 


m'' = 


(LA')i 72 


yuuv*^ / A' 0 0 

0 0 yctv*JA' 

0 ytcV*JA' yttVaJA' 

0 VdsV*^ /A' ydbVai /A' 

ysdV*JA' yssVaJA' 0 

0 0 ytbv* /A' 


(39) 


where 


UuiUj — UuiUj^iL/^ Cuj) t 

Vdidj = ydidjF{L,CQ,,Cdj) ■ 


(40) 


The up-type quark mass matrix is already block-diagonal. The diagonalization of the down-type mass matrix 
requires a more careful treatment. For the sake of simplicity, in the following analysis we denote the ij element of 
{mf') as m“ {mfA. The product of the down-type mass matrix and its adjoint 


^d 12 


73p 


771^2^22 






771^3 


1^21 p 


|W22P 


0 


777; 


d 12 
■331 


can be diagonalized in two steps: in first place, an approximate block diagonalization 


12 


^d ^d* 


~ I mf*mi 2 \m^ 

0 

is accomplished with the aid of the transformation matrix 


nd 12 


0 

0 



/ 

1 

0 

e 



0 

1 

0 


1 

-e* 

0 

1 


(41) 


(42) 


( 43 ) 
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and subsequently the diagonalization is completed through a unitary rotation of the upper block. This approximation 
is consistent provided Imfgl >> 1^121,1^131,1^221 |e| << 1. The resulting diagonalization matrices for the up 

and down sectors can be parameterized as 


Uu = 


Ud ^ 


/l 0 

0 

0 cos 9u 

sin due®"®®" 

\0 -sin6»„e-®‘^“ 

cos 9u 

/ cos 9d 

sin 9de'“f‘'‘ 

— sin dde"®®^"® 

cos 9d 

1 -e*cos6»d 

e* sindde®'®’'® 


(44) 


with 


tan 20U — ‘l\Z^\j 
tan 29 d = 2\Zd\/X~ 


(fiu = arg Zu , 
ipd = arg Zd , 


e = BdlAd , 


(45) 


and 


Xu = |w 
= Ito: 


u 12 
33l 
d 12 

22l 

12 


u |2 


1^32 r ± Im': 

Iwni P ± Ito 


u |2 


d 12 
12 l 5 


= 1 ^ 33 ! - - \m' 


12 1 




= TO^gTO^^ , 

Yd = , 

Bd = TO^gTO* . 




(46) 


Correspondingly, the quark mass eigenvalues can be expressed in terms of M^, defined in Eq. (l30l) . as 


TO,, = TO 


111 5 


md = ^M-{Xj,Yd), 

so that the CKM matrix is given by 


me = ^M-{X+,Yu), 
'ms = (Xj", Yd) , 


mt = j^XI+{X+,Yu) 


mb = Im^gl, 


VcKM = ulu, 

( 


(47) 


(48) 


cos 9d 


sin 9d e 

—cos 9u sin 6 d — sin 9u cos 9d€* cos 6 d cos 9u — sin 9u sin 9de* —6®“^“ sin 

cos 9u 


—e sindd sin^u — cos0„ cos^de* e ®‘^“ cos0dsin0„ — e®"^”® cos0„sindde* 


Hence, the quark sector Dirac CP phase (in PDG convention) and the Jarlskog invariant take the form 

^CP = arg(e) + (pd + ^u, 

Jqp ~ -|e| sin 2dd sin2du sin^cp . 


(49) 

(50) 


According to Eq. (HUl) . the size of up and down mass matrix elements is determined by the overlap of the 5D quark 
field zero mode profiles, i.e., to“ oc f^\L,CQ^)f^^\L,Cuj) and mf^ oc f)^\L,CQ.)f)^\L,Cdj)- If the wave function 
localization parameters cq. , , Cd, are chosen such that the quark zero mode profiles obey 


/f (L,ca) « f^°\L,cc) « /f (T,ct) 
(L, Cu) « /1j°^(T, ce) « (L, ct ), 

/®(L,Cd) « fld\LXs) « f^H\L,Cb) , 


(51) 


then the elements of to®® and to®^ approximately satisfy 


to“i < to“3 - TO^2 < ■^33 : ^12 ~ ^21 < "T.^2 < "^33 : 


Wi 3 < TO 33 , 


(52) 
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justifying the perturbative diagonalization performed on These relations imply that ^ \Yu,d\ holds, and 

therefore, a rough estimate for the mixing parameters and quark mass spectrum is 


^^23 

f^L\L,cc) 


m^2 


1^1 ^, 

mi3 

”133 

^ fi°\L,CT) ’ 


"122 

^ f''L\L,cc) ’ 

|e| ~ 

^33 


md 


ni I 


m.. 


rut ^ 
mb 




'■331 


m’ 


33 


Thus, in order to reproduce plausible quark masses and mixings, namely: 


(53) 


10 -1, 

Od- 

- 10-2, 

|e| 

- 10-3, 


mu : me 

: mt - 

- 10-1^ : 

: 10-2 : 

1 , 

(54) 

md'-ms \ 

: mb - 

- 10-3: 

:10-2 : 

1 , 



the quark zero mode profiles must observe the following hierarchy: 

/f (L,cc/) : f^j°\L,cc) : /f (^Tt) ~ IQ-^ : 10-' : 1, 

4°)(L,c,) : f^°\L,c,) : f^°\L,Ct) ^ IQ-^ : IQ-' : 1, 

Cd) : c,) : (T, Cb) ~ IQ-' : IQ-i : 1. (55) 


To conclude this section, an explicit realization of quark masses and mixings is presented. The choice cu = 1.97, 
cc = 1.92, CT = 1.83, Cu = -0.76, Cc = -0.62, c* = -0.56, Cd = -0.74, Cs = -0.69, Cb = -0.68, yuu = -0.438-0.954i, 
= -0.360 - 1.038i, ytc = 1.147 - 0.273i, yu = -0.372 - 1.073z, yds = -0.966 - 0.285i, ydb = 0.290 + 0.400i, 
ysd = 0.838 — 0.226i, ygg = —0.703 — 0.207i, ybb = 0.637 — 0.879i, generates the quark mass spectrum 


mu = 2.30 MeV, 
md = 4.80 MeV, 


me = 1.275 GeV, 
ms = 95.0 MeV, 


mt = 173 GeV, 
mb = 4.18 GeV, 


and fixes the magnitude of Vckm elements at 


I Vckm| 


/ 0.974 0.225 0.0035 \ 
0.225 0.973 0.0414 . 

^ 0.0089 0.041 0.999 / 


Finally, the obtained values for the Dirac GP phase and the Jarlskog invariant are 


(56) 


(57) 


(5^p = 1.25 , J^p = 3.06 X 10"^ . (58) 

The resulting quark masses and mixings are consistent with the current experimental data [l|, and the precision of 
the results can be improved by incorporating high order corrections, addressed in the next section. 


V. HIGH ORDER CORRECTIONS 

From the particle content and above transformation properties, one finds that nontrivial high order corrections 
to the charged lepton sector are absent in the present model. The next-to-leading order (NLO) corrections to the 
neutrino Yukawa interactions are given by 

= yG^[(rr)3'T2*¥,] i7vI/,^5(y-T) + h.c. (59) 

(A') 2 

However, the contribution of these terms to the neutrino masses and mixing parameters can be absorbed by a proper 
redefinition of the parameter j /22 after SSB. Hence, in order to estimate the effects of higher order corrections in this 
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sector, we need to investigate the Yukawa terms involving an additional (vyh/ h!)'^ suppression with respect to the 
lowest order terms in Eq. (HU, where we have introduced wir, to characterize the magnitude of ^ ~ . 

The contraction of the field products transforming as (3, -1,-1) under A(27) 0^4® Z4, with 

the flavon operators 




1 




(A')^ ' (A' 

as well as the combination of ~ (3,q—1) and 


(AO 


<<yl. 


(AO 


Trs‘^1 5 


1 


(AO^ 


1 


(A': 




(7 *2 


(AO 


(60) 


(61) 


provide the desired high order corrections to the neutrino Yukawa interactions. In the above expressions, the indices 
a, 6 = 0,1, 2 label the different singlets of A(27). Additional terms that can be absorbed into yn, 2/13, 2/22, 2/3i and 2/33 
have been omitted. Taking into consideration these corrections, the neutrino mass matrix rriv can be roughly written 
as 


rriv ~ 


(LAO^ V2 


A' 


2/11 A' 
0 

V 2/31 — 







2:21 

V 0 


0 X 23 

X 32 0 


(62) 


with Xij = XijF{L,ci,Ci,^)^ and Xij as dimensionless parameters of order 0(1). 

Working under the same numerical framework established in Section IlIIl one can readily estimate the shift in the 
neutrino oscillation parameters induced by high order corrections of the Yukawa interaction. Particularly, in Case I, 
taking Xij as random complex numbers with magnitudes ranging from 2 to 6, and uir = 0.1, the resulting deviations 
in the neutrino mixing parameters with respect to their LO values can be estimated as 


5s?2 ~ 0.01 Ssls - 0.01 <5s?3 ~ 0.001 SJcp 0.001. 


(63) 


On the other hand, the corrections to the neutrino mass splittings are negligible 


S{Amli) - 10“^eV^, S {A\ml^\) ~ 10”® eV^. 


(64) 


From Eq. ()63L it is clear that high order corrections can easily drive S42 into its Icr region while keeping the remaining 
parameters optimal. 

Turning to the quark sector, every bilinear formed by 4 /q^ and or 4/^^ can produce a high order correction 
to the Yukawa interaction whenever it is contracted with the adequate cubic flavon operator. Beside terms that can 
be absorbed by a redefinition of ymuj or ydidj in Eq. (1551) . all the NLO contributions can be classified into three 
categories: 

• Invariant products of Tt/T/Tc, ~ (lo,iA)—1) with 


I 


(CD 


1o,2^2 5 


1 


^1^2 • 


(65) 


Invariant products of Tyi/'I't, ^tH^s (lo,o,l)—1) with 


(A')^ 


■(?Olo,20’l ) 


7T70(er)i„4at, 

(AO 5 


(AO 




(AO 


^*3 


Invariant products of and ^ (lo,2,^—1) with 

-i- 1 -IrO 


(A' 


77T(Cr)lo4<^2% 

' I O 


(AO^ 


( 66 ) 


( 67 ) 
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Again, after symmetry breaking, the quark mass matrices rriu and rrid can be approximately written as 


m = 


m'^ = 




{LK')i 72 


yuuvl^jh! 0 0 

0 0 VcKJh! 

0 ytcvl^/J^' yttVo^jh! 

0 ydsvl^/h-' ydbVaJ^' 
yBdV*JA' yssVaJA' 0 

_V 0 0 ybbVaJA' 


/ ^IR 7 
V A' ) 


/ f IR Y 

V A' / 



( 68 ) 


Here we have defined = XuiUjF{L,CQ,,Cu^) and Xd^dj = A(L, cq,, c^Y, where the couplings Xu^uj and Xd^dj 
represent dimensionless parameters of order 0(1). As a numerical example, taking Xmujy Xd^dj as random complex 
numbers with magnitudes ranging from 1 to 4 for Xuc, Xcu, Xbd, Xcc, Xtu, Xdd, and from 2 to 6 for Xut, Xsb, Xbs, while 
keeping the values of CQ^, c„^, Cd^, yuiUj and ydid, reported in Section ITVl the order of deviation with respect to the 
LO values of the quark masses is 


(5to„ ~ 0.001 MeV, 
5md ~ 0.1 MeV, 


~ 10MeV, 
5ms ^ 0.1 MeV, 


5mt ~ 0.1 MeV, 
5mb ~ 0.5MeV. 


(69) 


The corresponding correction to the first order CKM matrix is of order 

/ 0.001 0.005 0.0001 
(5|Vckm| ~ 0.005 0.001 0.001 

\ 0.0005 0.001 0.00005 

and the values for the quark CP violating phase and the Jarlskog invariant are displaced by 

5{51p)^Q.1, <5J^p~10-7 


(70) 


(71) 


As for the lepton sector, it is not difficult to find parameter values reproducing the quark mass and mixing parameters 
required to fit the current experimentally observed values. 


VI. CONCLUSIONS 

We have proposed a hve-dimensional warped model in which all standard model fields propagate into the bulk. Its 
structure is summarized in the “cartoon” depicted in Figure [2] Mass hierarchies in principle arise from an adequate 
choice of the bulk shape parameters, while fermion mixing angles are constrained by relations which follow from 
the postulated A(27) flavor symmetry group, broken on the branes by a set of flavon fields. The neutrino mixing 
parameters and the Dirac CP violation phase are described in terms of just two independent parameters at leading 
order. This leads to stringent predictions for the lepton mixing matrix which should be tested in future neutrino 
oscillation experiments. Likewise the scheme also includes the quark sector, providing an adequate description of the 
quark mixing matrix. The effect of next-to-leading order contributions is estimated to be fully consistent with the 
experimental requirements. 
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FIG. 2: Pictorial description of the basic warped model structure, showing the UV (IR) peaked nature of the standard model 
fields. 


Appendix A: Group theory of A(27) and its representation 

The A(27) group is isomorphic to (Z3 (g) Z3) x Z3. It can be conveniently expressed in terms of three generators a, 
a' and b which satisfy the following relations: 

= a'^ = = 1, aa' = a'a, 

bab~^ = a~^a' ba'b~^=a. (Al) 

All A(27) elements can be written into the form with k,m,n = 0, 1, 2. The group has 11 conjugacy classes, 

given by 

IC^^) = {aa'^}, 

ICf ^ = {a^a'}, 

= {a, a', a^a'^}, 

dCg = {a , a , aa I, 

= {b,baa'^,ba^a'}, 

= {ba,ba',ba^a'^}, 

= {ba^,baa',ba'^}, 

acf’^^ = {62a,62a',&2a2^/2|^ 

acf= {b^a^, b^aa', 6^0'^} . (A 2 ) 


The A(27) has nine one dimensional representations, which we denote as lk,r {k, r = 0,1, 2), and two three dimen¬ 
sional irreducible representations 3 and 3. The explicit form of the group generators in each irreducible representation 
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Xio,o 

Xio,i 

Xlo,2 

Xii,o 

Xii,i 

Xll,2 

Xl2,0 

Xl2,l 

Xl2,2 

iC3 

Xs 

IGi 

1 

1 

1 

1 

1 

1 

1 

1 

1 

3 

3 

1 C)^> 

1 

1 

1 

1 

1 

1 

1 

1 

1 

36.2 

36 . 

IGf’ 

1 

1 

1 

1 

1 

1 

1 

1 

1 

36 . 

36.2 

3 ^ 

1 

OJ 


1 

UJ 

6.2 

1 

UJ 

6.2 

0 

0 


1 


u 

1 

6.2 

u 

1 

6.2 

UJ 

0 

0 


1 

1 

1 

CJ 

UJ 

UJ 

6.2 

6.2 

6.2 

0 

0 


1 

(Jj 


UJ 

6.2 

1 

6.2 

1 

UJ 

0 

0 


1 

a;2 

UJ 

UJ 

1 

6.2 

6.2 

UJ 

1 

0 

0 


1 

1 

1 


0.2 

6.2 

UJ 

UJ 

UJ 

0 

0 


1 

OJ 


UJ^ 

1 

UJ 

UJ 

UJ^ 

1 

0 

0 


1 


UJ 

0.2 

UJ 

1 

UJ 

1 

6.2 

0 

0 


TABLE VI: Character table of A(27). 


is 


Ifc.r ■ 

a = 6.’’, 


II 

e 


/ 6. 0 0 \ 


^ 6.2 0 0 \ 

3 : a = 

0 1 0 

, a! = 

0 6. 0 


^ 0 0 6.2 y 


1 0 0 1 / 


^ 6.2 0 0 ^ 


0 

0 

3 : a = 

0 1 0 

, a! = 

0 6.2 0 


0 

0 


U 0 \) 


b = 


b = 


b = 

/ 0 1 0 \ 

0 0 1 , 
Vl 0 0/ 

/ 0 1 0 \ 

0 0 1 , 

VI 0 o) 


(A3) 


where w = is the cube root of unity. Notice that 3 and 3 are complex representations dual to each other. From 

the character table of the group, shown in Table IVTl we can straightforwardly obtain the Kronecker products between 
the various representations 


0 — l[/c+fc'],[r+r'] 5 ^ ^ lfc,r — ^5 3 0 lfc,r — 3 , 

2 

3®3=^lfc,r, 3(g)3 = 3©3©3, 3©3 = 3©3©3, (A4) 

fc,r—0 

where [n] stands for n mod 3, whenever n is an integer. Starting from the representation matrices of the generators 
in different irreducible representations, we can calculate the Clebsch-Gordan (CG) coefficients for the Kronecker 
products listed above. All GG coefficients are presented in the form a ® fi, where stands for the elements of the 
first representation and /3j those of the second one. In the following, we adopt the convention q;[ 3 ] =00 = 0 ^ 3 - 

• l/j, j-/ = 


• 3 © = 3 


(ai) ©(,Si) =( aiPi ) 


/aiV 

02 

© 

(*). = 

^ a[i+r]/3i \ 

c'^ap+r]/?! 



^ ' lfc,r 

\ 6.2''a[3+^]/3i / 
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• 3 0 = 3 


^ ai'^ 

/ V 

<2[1—r]/3l \ 

<22 

0 ^4 = 

w''<2[2_r-]/3l 



Uj'^'^a[3_r]l3l J 


• 303 — X]fc,r=0 


<2l \ 



//3i\ 

<22 


0 

/32 

\a3 J 

3 


\^3 / 3 

(<2i/3i 

+ 

<22/32 

+ <23/33)10,0 

©(<23/3i 

+ 

<2i/32 

+ <22/33)10,1 

©(<22^ 

+ 

<23/32 

+ <21/33)10 2 


© ( q ! i / 3 i + 0j'^a2^2 + 0Ja3l33)x^ g 

© (q!3/3i + a;^ai/32 + a;a2/33)ii,i 
© (Q!2/3i + w^a3/32 + a;ai/33)ii 2 


© (ai/3i + a;a2^2 + w^q;3/33)i2^o 
© (a3/3i + wai/32 + a;^a2/33)i2,i 
© (a2/3i + wa3/32 + a;^ai/33)i2 ^ 


• 3 0 3 = 3si © 3 s2 © 3a 


^ <2i ^ 

\ 

// 3 i\ 

^ < 2 i/ 3 i \ ^ / 

< 22/32 ® 2 

<22 

0 

P 2 

\«3 y 

/ 

3 

\^3 / 3 

^ 3 / 33 / 3 ^^ V 


oisPi + ai/33 


3 So 


^ <22/33 ~ <23/32 
<23/3i — <2i/33 
\ ai/32 - <22/3 i J 


where the subscripts “5” and “o4” denote symmetric and anti-symmetric combinations respectively. 
• 3 0 3 = 3s, © 3 s2 © 3a 

^ ai 
<22 

\«3/3 



^ <2i/3i \ 
<22/32 

V ^3/33 j 3 


^ <22/33 + <23/32 \ 
<23/3i + ai /33 
\ <2i/32 + <22/3i / 

OSo 


^ <22/33 — <23/32 \ 
a3/3i - ai/33 
\ ai/32 - a2/3i / 3^ 


(AS) 


(A6) 


Appendix B: Vacuum Alignment 


In this Appendix, we shall investigate the problem of achieving the vacuum configuration in Eq. (1101) and Eq. (HID . 
For self-consistency, all fiavon fields </?, cti and 02 are treated as complex, given the form of the A(27) representation 
matrices, and the fact that the Z4 charge of <12 is purely imaginary. Since the flavons and cti, <T 2 are assumed 
to be localized at y = 0 and y = L respectively, the vacuum alignment problem is greatly simplified. At the UV 
brane ?/ = 0, the fiavon ip transforms in the manner listed in Table HI The scalar potential invariant under the flavor 
symmetry A(27) 0 Z4 0 Z4 can be written as: 


Vuv = M^{ifif*)i^,,+fi {<p<p )3 +/2 

LI- J lo,o L ^ -11 

+ /3 [(7’7’)3s, +/3* [(w)3s2 ( 35 * 35 *)3sJ 

L ^ -I In.o L ^ J In.o 


(Bl) 


with real couplings M^, fi and f 2 - Note that (p = {(pi,ip 2 ,ip 3 ) is a A(27) triplet 3, and its complex conjugate 
(p* = transforms consequently as 3. Focusing on the field configuration 


{‘P) = (1,1, 


(B2) 
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the minimum conditions for the UV potential read 
dV\jv i9Vuv dV\jv 


difl dip2 

leading to a non zero solution 


d^3 


+ 2(/i + f2 + fs + /a)! 


= 0 , 


\v^\ = - - 




(B3) 


(B4) 


2(/i + /2 + /a + /a) ’ 

that holds in a finite portion of parameter space with /i + /2 + /a + /a <0. 

Similarly, at the IR brane y = L, the most general renormalizable scalar potential biu involving the flavon fields 
CTi, (72 is 

+ [m3sA^*nss, 


+92 




93 


- lo,o 

+g4Cr?(7t^ +55Cr|cr^^ +56|crincr2p +57(?r)lo,okl|'‘ +58(Cr)lo,ok2| 

+59(er)lo,(^?+^?9*(er)lo,.'^r^ 

where all couplings, excluding and gg, are real. For this potential, the Case I alignment 


mssXnss, 


lo,o 

2 


■93 


mssAcnss, 


-I lo 


(O = (0,U5,0), (c7i)=- 


(ct2) = 


(B5) 


(B6) 


determines the minimization conditions 


c^hiR * * 2 r, 

= 0 , 

^Rir 


/'jl/r2 , o I |2 , I |2 , I |2\ 

— = + 2gi\v^\ +37^^11 +g3\Va2\ J 


dC. 

dViK 2 n 

= ggv,v^^=Q, 

dVm 
dal 
dVm 
dal 


= 0 , 


= Vcri(^M^^+2g4\vai\^ + g7\v^\^ + gelva^l^'^ =0, 




^^2 I + 1^2 


255 fa 


■ 58 fj 


■ ^6 Pcri 




= 0 . 


(B7) 


From the above equations, it is clear that non-trivial solutions in this sector are only achievable by fine tuning the 
gg parameter to satisfy gg = 0. This choice can be enforced by an additional dynamical mechanism capable of 
switching off the , cr? an d .,(7*^ terms in the potential. Such scenario could be naturally realized in a 


supersymmetric extension 20|, . As this possibility lies beyond the scope or the present work, we simply impose 


the condition gg = 0 in the general potential. Then, the obtained solutions are given by 

(g| - 4 g 4 g 5 )A 7 | + (25557 - gGgs)M^^ + (25453 - ggg 7 )Ml^ 

kd = — 


2(4515455 + 565758 - 5156 - 545! - 5557) 

2 _ (25557 - 5658)A7| + (5I - 45i55)M2^ + (25155 - 5758)A7, 


2(4515455 + 565758 - 5156 - 545! - 5557) 

2 _ (25458 - 5657)A7| + (25156 - g7gs)M^^ + (gf - 45i54)M, 


2(4515455 + 565758 - 5156 - 545! - 555?) 


(B8) 


The right-handed side of these expressions can be positive in a finite region of parameter space. Analogously, for the 
Case II vacuum configuration 


(0 = (1: l)f5 : (f^l) = f<Ti , ( 0 - 2 ) = fa 2 , 


(B9) 
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the minimization conditions are 

dViR 

dC2 

dViR 

'W 

dVm 
dat 
dVm 

dan 


= 

+ 2 (gi - 

b g 2 + W^ga -k wgg) ^ -b gr u 

|2 
(Ti 1 

1 |2 
+ ffS 1 ^0-2 1 

1 2 1 ’t' =('2 

+ ggVcT, + 

= 0 

= UIV^ 

M| + 2 (gi 

+ g2 + w^ga -b wgg) ^ + 

gr 

Vai 

2 1 
+ ^8P^72 

2 , 2 2 , 2 * *2 
+ ^ g9f(7i + ^ ggt^iTi 

= Vf 

+ 2 (gi - 

b g 2 + w^ga -b wgg) ^ -b g? u 

|2 
(Ti 1 

1 

+ ffs 1 ^0-2 1 

1 2 1 * =t'2 

+ ‘^ggv^i + ggv^^ 

= 0 

= Wo-1 

+ 2g4 

|2 1 |2 ^ 1 |2' 
^0-1 +56^^72 +3^7?;^ 

= 

0, 




= 

-^^2 + 2g5 

|2 1 |2 ^ 1 |2' 
+56^^71 +3^8^^ 

= 

0. 





= 0 , 


(BIO) 


Again, these equations are incompatible unless gg = 0. Once the coupling gg is enforced to vanish, we are left with 
three independent linear equations for the three unknown variables v^, Va-^ and The solutions can be easily found 
as 


= 


(4g4g5 - gi)Af| + (gggs - 2 g^gr)Ml^ + (ffgg? - 2 gigs)M^^ 
2g(gi - 4g4g5) + 6(g4gi + gsg? - gegrgs) 

(Sgegs - ^gb97)Ml + (dggs - 3gi)M^^ + (Sgygg - 2gge)Ml^ 

- 4g4g5) + 6(g4g8 + gsg? - gegrgs) 

(Sgeg? - 6g4g8)M^ + (Sgygs - 2gge)M^^ + (4gg4 - 
^gigi - 4g4g5) + 6(g4g| + gsg? - gegrgs) 


(Bll) 


where we have defined g = gi + g 2 + w^ga +wg|. Therefore, both (^) = (0,uj,0) and (^) = alignments 

can describe the local minimum of Vir, depending on the the parameter values. In the case of gi <C g, the VEV 
(0 = (0i j 0) is preferred over (^) = (1, w, l)uj, while Em is minimized by (^) = (1, w, l)uj for gi ^ g. 
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